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Amplitude and Phase of the Low-Frequency 
Ground Wave Near a Coastline 

James R. Wait 

A theoretical analysis is given for the amplitude and the phase change of the ground 
wave, originating from a distant transmitter on land, as it crosses a coastline. The land 
and sea are assumed to be smooth, and homogeneous with a sharp boundary of separation. 
Attention is focused on the effects that take place near the coastline when it is not permis- 
sible to employ arguments based on the principle of stationary phase. A limited comparison 
is made with the recent experimental work of Pressed', Ash well, and Fowler. 

1. Introduction 

With the increased use of low-frequency radio navigational systems, it has become im- 
portant to understand the mechanism of ground-wave propagation over inhomogeneously 
conducting terrain. Theoretical treatments of this problem usually employ an idealized two- 
part media having a sharp boundary separating homogeneous media. The extension of these 
analyses to three or more media is straightforward, although tedious. In most of the previous 
work, 1 the numerical values of the fields are only valid for distances from the boundary greater 
than a wavelength or so, as a consequence of the approximations employed. It is the purpose 
of the present note to investigate the field near the boundary by a refinement of the stationary 
phase method. 

2. Theoretical Development 

In an earlier paper (see footnote 1), analysis and computations were presented for (he 
propagation of ground waves over a mixed path on a flat earth. The boundary separating 
the two media of contrasting conductivity was assumed to be straight and approximately nor- 
mal to the propagation path. A relation was derived for the mutual impedance Z' between 
the terminals of the transmitting and the receiving antennas, both of which are represented 
by dipoles. The conductivity of the earth to the left of the boundary is denoted by o and that 
to the right by a lm It is convenient to let 

Z'=Z+AZ, (1) 

where Z is the mutual impedance for the two antennas, A and B, corresponding to an earth 
that is all homogeneous with conductivity c and AZ is the increment to account for the portion 
of the path to the right of the boundary with conductivity ci. From the early work it follows 
that 
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where /3=27r/free-space wavelength, and h a and h b are the effective heights of the antennas at 
A and B, At=47rX10 -7 , co=angular frequency, 7)=(i/jLoo/a) y2 . F(r +R , rj) is a slowly varying 
Sommerfeld attenuation function defined previously (see footnote 1), and is a function of the 
separation distance r +/^ and surface impedance rj. 



1 J. R. Wait, Mixed path ground wave propagation; 1. Short distances, J. Research NBS 57, 1 (1956) RP2687. References to the work of 
Millington, Clemmow, Bremmer, Furutsu, and others are given here. 
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Figure 1. Schematic representation of a land-sea boundary with a transmitter at A and a receiver at B. 



The corresponding expression for the increment AZ is an integration over the media to 
the right of the boundary, which, in figure 1, corresponds to positive values of x; 
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where r and R are the distances from the variable point P to A and B, respectively, and 8 is 
the angle between unit vectors at P in the direction of increasing r and R. The attenuation 
function F(r,r)) is of the same form as that in eq (2); however, F f (R,rj,rji) is some unknown 
slowly varying function that depends on R, rj, and ^i[=(i/xw/ <7 'i)^]- By utilizing a stationary 
phase principle, the area integration can be reduced to a line integral from x=0 to R , enabling 
a solution to be carried out. To this approximation, F f (R,y\ y r\i) could be replaced by F{R,rii). 
It was admitted that this procedure leads to a result that is not valid near the boundaiy. 
As a further consequence of this stationary phase argument, the reflected wave from the 
boundary is neglected. In the present instance, R is considered small compared with r . It 
then can be anticipated that the important contribution of the integration will be confined to 
a region where x and y are small compared with r . In this region it is permissible to assume : 

F(r, v )(*F(ro,n)^F(r +R , r,), 

F'iRtfjTix)^!, (subject to o-i >><?■) 

l/r^l/r , cos 8q±(x-R )/R, (3(r+R)^p(r +x+R). 

In other words, the boundary is in the Fraunhofer field of the source at A. Then it follows 
that 



AZ^x — 77) 



P 2 h a h b 
47r 2 r 



e-^oF(r 0j v) 



(4) 



238 



■JJ.'-E.^^O+i)** 



where 

R^[(x-R y+y 2 }x 
It is now noted tliat 
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which is a consequence of Campbell and Fosters 2 pair No. 917, which requires (3 to have a 
finite negative imaginary part. H ( g ) is the Hankel function of the second kind. Using this 
result, it follows that 
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where 



for JS ^>0, and 
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Integrating by parts and changing the integration variable leads to 

7=j/iI< 2 >(/3#o^ (9) 

for /?o>0, and 

I=%{HP(-pRo)-U^J~~ a U H^(z)e^dz\ (10) 

for 7?o<CO, where e is any arbitrarily small real number. 3 
Now it can be readily verified that 

~{ze^[HP(zmH^(z)]}=e^H^(2), (11) 

so the integrations can be carried to yield 

I=i J [(1 -ipRo)HI?>(pRo)-(iRoH? ! >(JiRo)] (12) 



/=* \ (l-if}RoW-pR o )+l3RoH™(-0R o )] (13) 



for //i,<0. After a little substitution and cancellation, it is easy to show that 

Z 1+ Z 

= l+<ie^e i *[{\-ia)Ht?\a)-aH'?\a)} (14) 



2 A. N. Campbell and R. Foster, Fourier integrals (D. Van Nostrand Co.. Inc., New York, N. Y., 1948). 

3 The finite value of c insures that the integrals converge in the Riemann sense. After^carrying out the evaluation, c can be set equal to zero. 
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with a=j&R >0, and 
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with a=fiR o y>0. In the above 
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where j mc is the frequency in megacycles per second; q is a measure of the conductivity contrast 
between the two media and is a positive or negative real number if displacement currents in 
the ground can be neglected. Also, as a consequence of the approximate boundary conditions, 
|?7K<?7o=1207r, and therefore |AZ/Z| is restricted to values small compared with unity. 

3. Discussion of Results 

It is instructive, at this stage, to discuss the behavior of the ratio Z' '/Z for large a (i. e., 
large electrical distances from the boundary). The asymptotic developments for the Hankel 
function can then be employed, to yield 

>+¥-+«<!)"[« Jri^ + •■ •] w> 

when R >0, and 

i+¥«^(^«~|>TB+"-] < 18 > 

when R <0. It is most significant to note that the right-hand side of eq (17), when the square- 
bracket term is replaced by unity, corresponds to the first term of a development in higher 
powers of V^CW^o) obtainable from the earlier paper (see appendix). That is, the formulas 
derived by the| stationary phase methods are probably only valid if a^>l corresponding to 
R large compared to the wavelength. On the other hand, the stationary phase development 
would yield a value of unity for the right-hand side of eq (18). This would also be valid at 
large distances, in wavelengths, from the boundary. The departure from unity can be inter- 
preted as a reflected wave from the boundary that forms, with the incident wave, a standing- 
wave pattern on the transmitter side of the boundary. A more detailed view of the situation 
is illustrated in figures 2 and 3, where the amplitude and phase of (1 + AZ/Z) are shown plotted 
as a function of R /X f° r various values of q. It is interesting to note that very near the bound- 
ary the field has a logarithmic singularity because of the behavior of the zero-order Hankel 
function as its argument approaches zero. It is believed that the solution is not strictly valid 
for very small values of a because the approximate boundary condition introduced in the earlier 
paper is not applicable at or close to a discontinuity in the surface impedance. 4 

The appropriate values of the amplitude and phase 1 + AZ/Z derivable by stationary phase 
methods are also shown in figures 2 and 3 and indicated by dotted curves. There are significant 
differences between the two methods of calculations for regions within a few wavelengths of 
the boundary. However, the phase curves on the seaside of the boundary are almost identical 
to stationary-phase evaluated case, except for regions very close to the boundary. 

It is interesting to compare these results with some phase measurements made very 
recently by Pressey, Ashwell, and Fowler. 5 They confirmed the presence of the phase recovery 
effect that is experienced by a wave passing from land to sea. They also observed systematic 
variations, about the median recovery effect, whose magnitudes are of the order of 5° very 



* The approximate boundary conditions are not valid for i?o/X less than about \vlvo\ or a small fraction of a wavelength for well-conducting 
ground. 

6 B. G. Pressey, G. E. Ashwell, and C. S. Fowler, Change of phase of a low frequency ground wave propagated across a coast line, Proc. 
Inst. Elec. Engrs. [B] 103, 527 (1956). 
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Figure 2. Relative amplitude of the field near the coastline. 
The ordinate is the relative amplitude of the field at B for a distant transmitter at A on the land (fig. 1). 
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Figure 3. Relative phase of the field near the coastline. 
The ordinate is the relative phase of the field at B (fig. 1), 
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near the coast and decay to a negligible value within a few wavelengths from the boundary 
toward the sea. These variations, which are not predicted by the present theory, are of a 
form that varied from path to path and seemed to be dependent on the topography of the 
coastline. In practically all cases, however, they observed a large and rather sudden increase 
of phase lag very near the boundary, which is in qualitative agreement with figure 3. There 
was also some evidence from Pressey, Ashwell, and Fowler's work that a standing wave does 
exist on the transmitter side of the boundary, although the data here is rather meager. 

4. Appendix 

The stationary phase approach gives (see footnote 1) 
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If r is large compared to R , it is permissible to use the approximation 

F(r +x 9 n) F(rp) 
(r +x)y> == r* 

and the power series expansion of the function F (R G —x, 771) to yield 



where 



F = Z =i Vi L u 2 3 + J 



Boulder, Colo., October 25, 1956. 
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